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We study the effect on CPP GMR of changing the order of the layers in a multilayer. Using a 
tight-binding simple cubic two band model (s-d) , magneto-transport properties are calculated in the 
zero-temperature, zero-bias limit, within the Landauer-Biittiker formalism. We demonstrate that 
for layers of different thicknesses formed from a single magnetic metal and multilayers formed from 
two magnetic metals, the GMR ratio and its dependence on disorder is sensitive to the order of the 
layers. This effect disappears in the limit of large disorder, where the results of the widely-used 
Boltzmann approach to transport are restored. 
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PACS: 73.23-b, 75.70-i, 75.70Pa 

Giant magnetoresistance (GMR) in transition metal 
magnetic multilayers is a spin filtering effect which 
arises when the magnetizations of adjacent layers switch 
from an anti-parallel (AP) to a parallel (P) alignment. 
The resistance in the anti-aligned state is typically higher 
than the resistance with parallel alignment, the difference 
being as large as 100%. This sensitive coupling between 
magnetism and transport allows the development of mag- 
netic field sensors with sensitivity far beyond that of con- 
ventional anisotropic magnetoresistance (AMR) devices. 
In the most common experimental setup, the current 
flows in the plane of the layers (CIP), and the resistance 
is measured with conventional multi-probe techniques. 
Measurements in which the current flows perpendicular 
to the planes (CPP) are more delicate because of the 
small resistances involved. Despite these difficulties the 
use of superconducting contacts || , sophisticated litho- 
graphic techniques || , and electrodeposition , makes 
such measurements possible (for recent reviews see refer- 
ences §§). 

A widely adopted theoretical approach to GMR is 
based on the semi-classical Boltzmann equation within 
the relaxation time approximation. This model has been 
developed by Valert and Fert |l(],[ll]], and has the great 
advantage that the same formalism describes both CIP 
and CPP experiments. In the limit that the spin diffu- 
sion length l s f is much larger than the layer thicknesses 
(ie in the infinite spin diffusion length limit), this model 



reduces to a classical two current resistor network, with 
additional possibly spin-dependent scattering at the in- 
terfaces [ p2[ . Despite the undoubted success of this de- 
scription recent experiments Jl3|,[l4| have drawn atten- 
tion to the possibility of new features which lie outside 
the theory. Two important and central predictions of this 
model are that the CPP GMR ratio is independent of the 
number of bilayers in the case that the total multilayer 
length is not constrained to be constant, and furthermore 
is independent of the order of the magnetic layers in the 
case of different magnetic species. An apparent viola- 
tion of the first prediction has been observed in CIP and 
CPP measurements fll^JlSfl, a nd of the second prediction 
in CPP measurements |l3|]l4f | . However a convincing the- 
oretical explanation is lacking. 

The aim of this letter is to provide a quantitative 
description of the breakdown of the resistor model in 
diffusive CPP multilayers in the limit of infinite spin- 
relaxation length. To illustrate this breakdown, con- 
sider a multilayer consisting of two independent building- 
blocks, namely a (N/M) and a (N/M') bilayer, where M 
and M' represent magnetic layers of different materials or 
of the same material but with different thicknesses and N 
represents normal metal 'spacer' layers. From an experi- 
mental point of view M and M' must possess different co- 
ercive fields, in order to allow AP alignment. In the case 
of |l4|] this is achieved by considering respectively Co and 
Ni84Fei6 layers with Ag as non-magnetic spacer, while in 
p3[ both the layers are Co (with Cu as spacer) but with 
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different thicknesses (respectively lnm and 6nm). Two 
kinds of multilayer can be deposited. The first, that we 
call type I, consists of a (N/M/N/M')x/z sequence where 
the species M and M' are separated by an N layer and 
the group of four layers is repeated fi times. The second, 
that we call type II, consists of a (N/M)x/i(N/M')x/i se- 
quence, where the multilayers (N/M)x/i and (N/M')x/i 
are arranged in series. If the coercive fields of M (Hm) 
and M' {Hm 1 ) are different (eg Hm < Hm 1 ) and if N 
is long enough to decouple adjacent magnetic layers, the 
AP configuration can be achieved in both type I and 
type II multilayers by applying a magnetic field H whose 
intensity is Hm < H < Hm 1 - The AP configuration 
is topologically different in the two cases, because in 
type I multilayers it consists of AP alignment of adja- 
cent magnetic layers (conventional AP alignment), while 
in type II multilayers it consists of the AP alignment be- 
tween the (N/M)x/x and (N/M')x/x portions of the mul- 
tilayer, within which the alignment is parallel (see figure 
|l|a and figure 0b). From the point of view of a resistor 
network description of transport, the two configurations 
are equivalent, because they possess the same number of 
magnetic and non-magnetic layers, and the same num- 
ber of N/M and N/M' interfaces. Hence the GMR ratio 
must be the same. In contrast the GMR ratio of type I 
multilayers is found experimentally to be larger than that 
of type II multilayers , and the difference between 

the two GMR ratios increases with the number of bilay- 
ers. In the case of jl3| the GMR ratio of both type I and 
type II multilayers increases with the number of bilayers, 
which again lies outside the resistor network model. 

In this Letter we demonstrate for the first time that 
a description which incorporates phase-coherent trans- 
port over long length scales can account for such exper- 
iments. To illustrate this we have simulated type I and 
type II multilayers using a Co/Cu system with differ- 
ent thicknesses for the Co layers, namely icu = 10AP, 
< C o = 10AP, t' Co = 40 AP. The technique for comput- 
ing transport properties is based on a three dimensional 
simple cubic tight-binding model with nearest neighbor 
couplings and two degrees of freedom per atomic site. 
The general spin-dependent Hamiltonian is 

where a and (3 label the two orbitals (which for conve- 
nience we call s and d), i,j denote the atomic sites and a 
the spin, ef 17 is the on-site energy which can be written as 
ef = Eq + uh8 a d with h the exchange energy and a = — 1 
((7 = +1) for majority (minority) spins. In equation (Q), 
l^f a = is the hopping between the orbitals a and 

(3 at sites i and j, and <f ai (c^\) is the annihilation (cre- 
ation) operator for an electron at the atomic site i in an 
orbital a with a spin a. h vanishes in the non-magnetic 
metal, and 7"^ is zero if i and j do not correspond to 



nearest neighbor sites. Hybridization between the s and 
d orbitals is taken into account by the non- vanishing term 
7 sd . We have chosen to consider two orbitals per site in 
order to give an appropriate description of the density of 
states of transition metals and to take into account inter- 
band scattering occurring at interfaces between different 
materials. The DOS of a transition metal consists of a 
narrow band (mainly <i-like) embedded in a broader band 
(mainly sp-like). This feature can be reproduced in the 
above two band model, as shown in reference |l6| ] , where 
the appropriate choices for 7™^ and ef in Cu and Co are 
discussed. 

We analyze the simplest generic model of disorder, in- 
troduced by Anderson within the framework of the lo- 
calization theory jlTj , which consists of adding a random 
potential Vi to each on-site energy, with a uniform dis- 
tribution of width W (-W/2 < V < W/2), centered on 
V = 

?r = <?" + V- (2) 

The conductances and GMR ratios are calculated 
within the Landauer-Buttiker theory of transport |Ig| ] 
using a technique already presented elsewhere fL9[ . In 
figure H we present the mean GMR ratio for type I (type 
II) multilayers GMRi (GMRn) and the difference be- 
tween the GMR ratios of type I and type II multilay- 
ers AGMR=GMRi-GMR n , as a function of n for differ- 
ent values of the on-site random potential. The average 
has been taken over 10 different random configurations 
except for very strong disorder where we have consid- 
ered 60 random configurations. In the figure we display 
the standard deviation of the mean only for AGMR be- 
cause for GMRi and GMRn it is negligible on the scale 
of the symbols. It is clear that type I multilayers pos- 
sess a larger GMR ratio than type II multilayers, and 
that both the GMR ratios and their difference increase 
for largeu. These features are in agreement with exper- 
iments jn|[n| and cannot be explained within the stan- 
dard Boltzmann description of transport. The increase of 
the GMR ratio as a function of the number of bilayers is 
a consequence of enhancement of the spin asymmetry of 
the current due to disorder. In fact, despite the Anderson 
potential being spin-independent it will be more effective 
on the d band than on the s band, because the former 
possesses a smaller bandwidth. Since the minority spin 
sub-band is dominated by the d-electrons and the ma- 
jority by the s-electrons, the disorder will suppress the 
conductance more strongly in the minority band than 
in the majority. Moreover, since the transport is phase- 
coherent, the asymmetry builds up with the length, re- 
sulting in a length-dependent increase of the GMR ratio. 
The different GMR ratios of type I and type II multi- 
layers can be understood by considering the inter-band 
scattering. Both multilayers possess the same conduc- 
tance in the P alignment, while the conductance of type 
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I multilayer in the AP alignment is smaller than that of 
type II. The inter-band scattering is very strong when an 
electron crosses phase-coherently a region where the mag- 
netizations have opposite orientations, and this occurs in 
each (N/M/N/M') cell for type I multilayers, while only 
in the central cell for type II multilayer (see figure [l]a and 
|l|b). Hence the contribution to the resistance in the AP 
alignment due to inter-band scattering is larger in type 

I than in type II multilayers. Finally when the elastic 
mean free path is comparable with a single Co/Cu cell 
one expects the resistor model to become valid. To illus- 
trate this feature, figure || shows that in the case of very 
large disorder (W = 1.5eV), AGMR vanishes within a 
standard deviation as predicted by the Valert and Fert 
theory. 

As a second example in which the dependence of the 
GMR ratio on disorder changes when the multilayer ge- 
ometry is varied, consider the system whose AP align- 
ment is sketched in figure 0c and [l]d. In this case M and 
M' are different materials chosen in such a way that the 
minority (majority) band of M possesses a good align- 
ment with the majority (minority) band of M'. Moreover 
the thickness of the N layers has been chosen in order to 
allow an AP alignment of the magnetizations of adjacent 
magnetic layers in both type I and type II multilayers. In 
this case both type I and type II multilayers exhibit con- 
ventional P and AP alignments, but their potential pro- 
file is quite different. In figure ^ we present a schematic 
view of the potential profiles for type I and type II multi- 
layers for both the spins in the P and AP configuration. 
A high barrier corresponds to large scattering and a small 
barrier corresponds to weak scattering. The dashed line 
represents the effective potential for material M and and 
the continuous line for material M'. Figure [| illustrates 
that type I multilayers possess a high transmission spin- 
channel in the AP alignment, and hence the resulting 
GMR ratio will be negative. In contrast type II multi- 
layers do not possess a high transmission channel (there 
are large barriers for all spins in both the P and AP con- 
figuration) and the sign of the GMR ratio will depend 
on details of the band structure of M and M'. Consider 
the effects of disorder on these two kinds of multilay- 
ers. Using the same heuristic arguments as above we 
expect that the GMR ratio of type I multilayers will in- 
crease (become more negative) as disorder increases, in 
the case of disorder that changes the spin asymmetry of 
the current. This is a consequence of the fact that, in 
common with the conventional single-magnetic element, 
one of the spin sub-bands in the AP alignment is domi- 
nated by weak s-electrons (small barrier), which are only 
weakly affected by disorder. It is clear that this system 
is entirely equivalent to conventional single-magnetic el- 
ement multilayers discussed above. In contrast for type 

II multilayers there are no spin sub-bands entirely domi- 
nated by the weak scattering (small barriers) s-electrons, 
and all spins in either the P and AP configuration will 



undergo scattering by the same number of high barriers. 
In this case the effect of disorder will be to increase all the 
resistances and this will result in a suppression of GMR. 
Moreover it is important to note that in the completely 
diffusive regime, where the resistances of the different 
materials may be added in scries, the GMR ratio will van- 
ish if ~ R^, , where R£ is the spin-dependent 

resistance of the material A. To verify this prediction 
we have simulated both type I and type II multilayers 
using the parameters corresponding to Co and Fe72V28 
of reference respectively for M and M', and cor- 

responding to Cu for N. This choice was motivated by 
the fact that a reverse CPP-GMR has been obtained for 
(Fe 7 2V 28 /Cu/Co/Cu)x A i multilayers @. The GMR ra- 
tio for type I and type II multilayers is shown in figure |^, 
which illustrates the remarkable result that the GMR ra- 
tio of type I multilayers increases with disorder, while for 
type II structures it decreases. As explained above this 
is due to an enhanced asymmetry between the conduc- 
tances in the P and AP alignment for type I multilayers, 
and to a global increase of all the resistances for type II 
multilayers. As far as we know there are no experimental 
studies of the consequences of the geometry-dependent 
effect described above, and further investigation will be 
of interest, in order to clarify the role of the disorder in 
magnetic multilayers. 

Despite the fact that GMR was discovered more than 
ten years ago, it continues to present fascinating insights 
into transport in magnetic heterostructures. In this Let- 
ter we have addressed a new issue which lies outside the 
widely-adopted Boltzmann description of GMR, namely 
that changing the order of magnetic multilayers can sig- 
nificantly alter the magnetoresistance [ p~3|Jl4| . We have 
shown that this effect is a consequence of phase coherence 
on a length scale greater than the layer thicknesses. 
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FIG. 1. AP configuration for type I and type II multilayers 
described in the text. Figures (a) and (b) refer to the cases 
of thick and thin Co layers with different coercive fields. Fig- 
ures (c) and (d) refer to different magnetic metals coupled 
through non-magnetic spacers via exchange coupling. The 
black blocks represent Cu, the white Co and the hatched 
Fe72V28- The arrows indicate the direction of the magne- 
tizations. Note the difference in the case of figure (b) where 
the AP alignment occurs between the two halves of the mul- 
tilayer. 
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FIG. 2. GMR for type I (a) and type II (b) multilayers, 
and AGMR (c) in the case of thin (10AP) and thick (40AP) 
Co layers, as a function of the number of double bilayers 
Co/Cu/Co/Cu for different values of disorder. The symbols 
represent respectively W = (•), W = 0.3eV (□), W = 0.6eV 
(x), W = 1.5eV (O). 
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FIG. 3. Heuristic scattering profiles for type I (a) and type 
II (b) multilayers of the second example discussed in the text. 
The dashed and continuous lines represent respectively scat- 
tering potentials of material M and M'. 
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FIG. 4. Different geometry-induced behavior of the GMR 
ratio as a function of disorder in multilayers composed of Co 
and Fe72V28- In this case all the layer thicknesses are fixed 
at 10AP. The open (closed) symbols represent type I (type 
II) multilayers discussed in the text. The circles are the dis- 
order free case, squares and diamonds are for random on-site 
potentials of 0.6eV and 1.2eV respectively. 
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